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ABSTRACT 

We  introduce  a  new  type  of  paxtition  called  a  parallel  planes  partition.  We  prove  there  exists 
a  parallel  planes  partition  of  any  set  of  n  points  in  arbitrary  dimension.  This  partition  yields 
a  data  structure  for  the  half-space  retrieval  problem  in  arbitrary  dimension;  it  has  linear  size 
and  achieves  a  sublinear  query  time.  Also,  we  give  efficient  algorithms  for  computing  this 
partition. 


This  work  was  supported  in  part  by  NSF  grant  DCR-84-01633  and  by  an  IBM 
faculty  development  award. 


1.  Introduction 

The  half-space  retrieval  problem  is  the  following.  Given  a  set  of  n  points  in  d- 
dimensional  Euclidean  space,  preprocess  them  so  as  to  be  able  to  quickly  answer  the  query; 
how  many  points  lie  in  the  query  half  space  H.  (A  variant  of  the  problem,  the  listing  prob- 
lem, is  to  ask  for  a  list  of  the  points  in  the  half  space.)  It  is  assumed  that  many  such  queries 
will  be  made.  Thus  it  is  reasonable  to  preprocess  the  set  of  points  and  to  amortize  the  cost  of 
diis  preprocessing  over  the  (many)  queries.  Our  Hrst  concern  here  is  with  the  query  time 
and  the  space  used  by  the  data  structiue  for  holding  the  preprocessed  information.  We  note 
that  the  naive  seardiing  algorithm  takes  linear  time.  Thus,  we  aim  for  a  sublinear  search 
time  and  a  data  structure  using  linear  space.  A  second  concern  is  to  find  relatively  efficient 
algorithms  for  building  the  data  structures. 

Recently,  an  elegant  approach  to  this  problem  was  discovered  by  Willard  [W].  There 
are  two  basic  lemmas  underlying  his  construction.  Hrst,  given  a  set  of  n  points,  the  plane 
can  be  partitioned  by  two  straight  lines  into  4  quadrants,  so  that  each  (open)  quadrant  holds 
at  most  r/4  points.  This  leads  to  a  recursive  storage  of  the  points  in  a  4-way  balanced  tree 
(ignoring,  for  now,  points  lying  on  the  partitioning  lines).  Second,  any  line  intersects  only  3 
of  the  4  quadrants  defined  by  the  two  partitioning  lines.  This  implies  that  in  carrying  out  a 
half-space  query  at  most  three  of  the  four  subtrees  of  the  root  need  be  explored  further. 
Recursive  application  of  this  observation  yields  a  sublinear  search  time. 

By  tuning  his  data  structure  (dividing  the  set  into  6  equal  parts  with  3  lines,  rather  than 
4  parts  with  2  lines)  Willard  obtained  a  linear-sized  data  structure  that  supported  a  query 
time  of  0{n'^'').  (Henceforth  we  wiU  not  refer  to  the  size  of  the  data  structure,  for  it  is 
always  linear  in  this  paper.)  Subsequent  work  by  Edelsbrunner  and  Welzl  [EW]  improved 
the  query  time  to  0{n'^^^). 

Using  the  same  approach  Yao  [Y1,Y2]  obtained  a  similar  result  in  3  dimensions.  She 
proved  that,  given  any  set  of  n  points,  the  space  can  be  divided  by  3  planes  into  8  regions,  so 
that  each  open  region  contains  at  most  one  eighth  of  the  points.  This  yielded  an  algorithm 
for  the  half -space  retrieval  problem  with  a  sublinear  query  time.  Dobkin  and  Edelsbrunner 
[DE1,E]  then  further  improved  the  algorithm  that  is  derived  from  this  partition,  obtaining  a 
query  time  of  0(n°  *'). 

The  natural  generalization  in  d  dimensions  is  to  seek  to  partition  the  space  into  2'' 
regions,  so  that  eadi  open  region  contains  at  most  n/2''  points.  Such  a  partition  would  yield 
an  algorithm  for  half  space  retrieval,  with  sublinear  query  time,  similar  to  Willard's.  It  is 
natural  to  try  to  achieve  this  partition  with  d  planes.    However,  Avis  [A]  showed  that  for 


d^S  there  exist  sets  of  n  points  for  which  there  is  no  such  partition.  In  fact  he  showed  that, 
in  general,  any  partition  by  planes  requires  at  least  {2''-l)/d  planes.  While  not  proving  that 
there  is  no  algorithm  with  a  sublinear  search  time,  this  result  does  have  a  somewhat  negative 
flavor. 

Recently  A.  Yao  and  F.  Yao  [YY]  showed  that  there  exists  a  partition  in  d  dimensions, 
for  any  rf.  Their  partition  divides  the  space  into  2**  regions,  using  parts  of  S^**"^-!  planes, 
so  that  each  of  the  open  regions  in  the  partition  contains  at  most  L^**  of  the  points.  We  show 
a  similar  resiilt  for  the  parallel  planes  partition:  it  divides  the  space  into  2''  regions,  using 
parts  of  2**"^  planes,  so  that  each  of  the  open  regions  in  the  partition  contains  at  most  l;^''  of 
tije  points.  This  is  a  generalization  of  the  partition  described  in  [CI].  It  is  not  a  variant  of 
the  partition  given  in  [YY];  the  underlying  ideas  are  similar,  however. 

Both  these  partitions  immediately  yield  a  linear  sized  data  structure  for  the  half-space 
retrieval  problem  in  d  dimensions,  supporting  a  sublinear  query  time.  (The  parallel  planes 
partition  yields  a  slightly  worse  query  time  than  the  partition  in  [YY].)  Further  applications 
of  the  partition  result  include  the  circle  retrieval  problem  [Yl]  and  other  query  problems 
[DE2].  Our  work  and  [YY]  extend  these  results  to  arbitrary  dimension;  for  example,  both 
partitions  provide  a  data  structure  for  the  sphere  retrieval  problem  supporting  a  sublinear 
query  time,  in  arbitrary  dimension. 

An  interesting  aspect  of  our  work  is  that  the  partition  is  not  by  d  planes  (in  fact  it  uses 
parts  of  2'*"^  planes).  So  the  resuJt  proved  by  Avis  does  not  rule  out  our  construction.  TTus 
opens  up  the  prospect  of  finding  other  partitions  for  point  sets  in  arbitrary  dimension.  In 
fact,  it  suggests  one  might  seek  other  partitions  enjoying  the  following  properties: 

(a)  Every  open  region  defined  by  the  partition  contains  at  least  some  constant  fraction  of 
the  points.  (Or,  more  precisely,  the  open  region  plus  the  separating  surfaces  contain  at 
least  some  constant  fraction  of  the  points.) 

(b)  For  any  plane,  there  is  at  least  one  region,  deHned  by  the  partition,  that  the  plane  does 
not  intersect. 

A  partition  satisfying  property  (b)  is  called  an  avoiding  partition  (any  plane  avoids  at  least 
one  region).  A  partition  satisfying  both  properties  is  called  a  separating  partition  (property 
(a)  was  first  proposed  in  [Yl]).  A  separating  partition  immediately  yields  an  algorithm  for 
half  space  retrieval,  of  the  form  given  by  Willard,  having  a  sublinear  query  time.  A  caveat 
should  be  made  here.  In  Willard's  data  structure  it  is  necessary  to  be  able  to  perform  half 
space  queries  on  the  separating  lines,  if  they  contain  any  of  the  points.  Similarly,  to  obtain 
an  efficient  algorithm  from  a  class  of  separating  partitions,  we  need  to  be  able  to  perform 


half  space  queries  on  the  separating  surfaces,  and  these  queries  must  be  performed  faster 
than  the  general  queries.  (Our  partition  uses  (^— l)-dimensional  planes  for  the  separating 
surfaces,  and  for  these  surfaces  the  requirement  is  met.) 

In  the  second  part  of  the  paper  we  give  algorithms  for  computing  the  parallel  planes 
partition.  The  algorithms  compute  a  partition  in  d^3  dimensions  in  time  Oin\o^^''^n),  where 
f(d)  =  2'^~^—d.  (The  previous  best  result  computed  a  weaker  partition  in  da 3  dimensions 
in  time  0{n'')  [YY];  in  2  dimensions,  Megiddo  [M2]  gave  a  linear  time  algorithm.)  Our  algo- 
rithms are  based  on  Megiddo's  ingenious  technique  for  using  a  parallel  algorithm  for  one 
problem  to  derive  an  efficient  serial  algorithm  for  a  related  problem  [Ml].  In  our  case  the 
parallel  algorithm  will  be  for  Ae  sorting  problem,  so  we  are  able  to  use  the  enhancement  of 
Megiddo's  technique  given  by  Cole  [C2].  Our  algorithms  will  require  repeated  recursive  use 
of  Megiddo's  technique;  this  was  illustrated  previously  in  [CSY]. 

These  algorithmic  techniques  can  also  be  used  to  obtain  algorithms,  with  the  same  run- 
ning times,  for  computing  a  special  case  of  the  partitions  described  in  [YY];  this  partition 
uses  (parts  of)  2*^'^  separating  planes. 

In  section  2  we  describe  the  parallel  planes  partition  and  prove  that  it  is  avoiding.  In 
section  3  we  show  that  it  is  separating.  In  section  4,  we  describe  several  algorithms  for  com- 
puting the  parallel  planes  partition. 

2.  The  PartltloD  Is  Avoiding 

We  make  the  notion  of  a  partition  precise. 

Definition:  [Yl]  Let  5  be  a  set  of  n  points  in  </-dimensional  space.  A  partition  11  for  ,S  is  a 
pair  (R,P),  where  R  is  a  set  of  disjoint  open  regions  and  P  is  a  set  of  (J- l)-dimensional  sur- 
faces, such  that  every  point  of  5  either  lies  in  a  region  /?  c  R,  or  on  a  surface  />  e  P.  A 
region  /{  e  R  is  said  to  be  a  hn  — region  with  respect  to  II  if  at  least  8n  points  of  S  are  con- 
tained in  /?  y  P  (i.e.  the  remaining  open  regions  of  R  contain  at  most  (1-8)1  points).   We 

say  n  =  (R,P)  is  a  hn- partition  if  every  /?  e  R  is  a  8n-region. 

We  describe  the  parallel  planes  partition  next.  As  an  aid  to  our  intuition,  we  start  with 
the  partition  11  in  3  dimensional  space.  We  use  4  partitioning  planes:  P,  Q,  R^^,  R^-  Ri  and 
R2  are  parallel;  further,  we  only  use  parts  of  the  planes  R^  and  R2,  as  illustrated  in  figure  1 
(hereafter,  when  we  refer  to  a  plane,  we  intend  that  part  of  the  plane  in  the  partition).  It  is 
clear  that  these  planes  defme  8  disjoint,  open  regions. 


Lemma  1:  Any  plane,  M,  intersects  at  most  seven  of  the  regions  defined  by  n. 

Proof:  The  approach  used  in  this  proof  was  suggested  by  Edelsbrunner  [E].  We  show  that 
the  planes  P,  Q,  R^,  J?2,  define  at  most  seven  faces  on  M.  We  can  then  conclude  that  M 
intersects  at  most  seven  of  die  regions  defined  by  II. 

P  and  Q  defme  four  open  regions  in  M;  we  refer  to  these  regions  as  quadrants.  Let  R 
be  parallel  to  R^  (and  hence  to  R2)  and  suppose  that  R  passes  through  PnQHM  (see  figure 
2).  Then  R  does  not  intersect  at  least  one  pair  of  opposite  quadrants.  Only  one  of  the  planes 
Ri  and  Rj  can  intersect  either  of  the  quadrants  in  this  pair;  without  loss  of  generality  suppose 
that  R^h  this  plane.  But  R^  can  intersect  at  most  one  of  the  quadrants  in  this  pair,  since  R^  is 
parallel  to  R.  We  deduce  that  P,  C.  ^1.  and  R2  define  at  most  seven  faces  on  Af .  o 

In  d  dimensions,  the  partition  II  uses  2'''^  planes  (or  rather,  parts  of  these  planes).  We 
use  planes  Pj,  im^d,  where  for  i>l,  0s;<2'~^,  and  for  j  =  1,  j  =  0.  We  say  a  plane  Pj  is 
at  Uvel  i.  The  planes  are  each  given  an  orientation.  We  label  the  open  regions  defined  by 
the  planes  at  levels  1  through  i,  as  follows.  Uj  is  such  a  region;  the  k-th  digit  of  7,  in  the 
binary  representation  (counting  from  the  right  end),  is  equal  to  1  if  the  region  is  on  the  posi- 
tive side  of  its  bounding  level  *  plane,  and  to  0  if  it  is  on  the  negative  side.  We  note  that  l/j 
is  partitioned  by  a  level  i+1  plane  into  l/^j  and  L^*^.  The  signature  of  region  C/j  is  defined 
to  be  the  i-bit  number  j.  Also,  the  signature  of  any  point  in  Uj~^  is  defined  to  be  the 
(d-  l)-bit  number  j. 

The  partition  II  in  3  dimensions  was  described  above.  The  partition  11  ind  dimensions, 
d>3  is  obtained  as  follows.  Let  jcj  be  an  axis.  We  require  the  level  t  planes,  i<d,  to  be 
parallel  to  x^.  We  also  require  the  level  1  planes,  i<d,  to  form  a  (<f-l)-dimensional  parti- 
tion in  any  (J- l)-dimensional  subspace  perpendicular  to  x^.  It  remains  to  define  the  levd  d 
planes.  Plane  Pj  is  chosen  to  divide  the  opposite  regions  L^~^  and  U^~^,  where  a  =  2j,  and 
P  =  a,  the  complement  of  the  (J-l)-bit  ntmiber  a;  the  parts  of  the  plane  not  intersecting 
these  two  regions  are  discarded.  We  will  also  impose  certain  parallelism  constraints  on  these 
planes;  we  describe  these  below. 

A  remark  on  notation:  We  abuse  notation  by  using P'^  f)  P*  f]  ■••  f)  P^  to  denote  an  arbi- 
trary plane  parallel  to  some  P/  Cl  P/  Ci  '"  Ci  ^}>  wherever  used,  it  wiU  be  well  defined, 
as  shown  later. 

Definition:  Let  t/  be  a  region  and  let  P  and  Q  be  bounding  planes  for  t/.  U  is  mixed  with 
respect  to  P  and  Q  if  it  is  on  the  positive  side  of  one  of  these  planes  and  on  the  negative  side 


of  the  other;  it  is  unmixed  if  (/  is  on  the  same  side  of  both  of  these  planes. 
Fact:   For  each  plane  Pj,  3s/srf,  0sj<2'~^  the  two  regions  divided  by  Pj  have  the  same 
bounding  planes;  also,  they  are  either  both  mixed  or  both  unmixed,  with  respect  to  any  pair 
of  bounding  planes,  one  plane  at  level  k,  the  other  at  level  l^k,  l,k<  i. 

The  parallelism  is  enforced  in  stages. 

For  the  0-th  stage,  all  the  level  d  planes  are  constrained  to  be  parallel  in 
^^  n  P^  n  •••  n  ^''"^-  We  denote  the  set  of  level  d  planes  by  5g. 

The  i-th  stage,  l^isd-3,  proceeds  as  follows.  Suppose,  at  the  end  of  the  (j-l)-th 
stage,  we  have  obtained  sets  5/"'"^^  0^j<2'~K  We  separate  the  planes  in  S/~'*^  into  two 
sets:  Sff'  and  Sff^^.  A  plane  is  placed  in  S^~'  (resp.  5^"+i)  if  it  divides  unmixed  (resp. 
mixed)  regions  U''~^,  where  we  consider  the  regions  with  respect  to  their  bounding  planes  at 
levels  d—i-1  and  d-i.  The  planes  in  Sff'  are  constrained  to  be  parallel  in 
P^  O  P^  0  ••   n  /*'*~'~^  as  are  the  planes  in  5^"+'i. 

We  also  have  a  vacuous  (J-2)-th  stage,  solely  to  define  the  sets  Sj,  0^j<2''~^;  each  of 
these  sets  contains  just  one  plane. 

We  can  compute  similar  sets  for  the  level  k  planes,  3^k<d;  we  denote  these  sets  by 
S^"'-*  (and  extend  the  notation  to  Jk  =  i).  We  note  that  d-i^k.  We  remark  that  the  planes 
in  5y ''•*  are  exactly  the  level  k  bounding  planes  for  the  regions  divided  by  the  planes  in 
^27+t*''  O^^^^"",  l^r^d-k.  (This  is  easy  to  check,  by  induction  on  i,  for  the  division  of 
^2'7+;*'  into  two  sets  is  made  according  to  whether  the  planes  in  Sj"^-**'  divide  mixed  or 
unmixed  regions  with  respect  to  the  bounding  planes  at  levels  d-i  and  d-i+1.) 
Lcmnui  2:  P^  f^  •••  (^  P'  ii  well  defmed.   (And  hence  the  partition  II  is  well  defined.) 

Proof:  The  result  follows  by  a  simple  induction  on  /.  It  suffices  to  note  that  stage  0  of  the 
construction  above  ensures  the  result.  Q 

Remark:  (This  remark  is  relevant  to  section  3.)  We  define  the  signature  of  set  Sf''  as  fol- 
lows. Osnsider  the  regions  divided  by  planes  in  Sj~'.  Let  U''~^  be  the  region  on  the  negative 
side  of  its  level  j  bounding  planes,  for  l^j^d-i-l  (it  is  easy  to  show  that  there  is  exactly 
one  such  region).  The  number  comprising  the  first  i  bits  (from  the  left  end)  of  the  signature 
of  l/*'^  is  defined  to  be  the  signature  of  Sf~'.  This  signature  is  denoted  <t(J). 

We  now  show  that  the  partition  is  avoiding.    Some  notation  will  be  helpful.   Let  M  be 
an  arbitrary  plane.   Let  P  be  a  level  d  plane.   Let  A  and  B  be  the  two  regions  divided  by  P. 


Note  that  A  and  B  have  the  same  bounding  planes.  Define  Op  to  be  the  point  that  is  the 
intersection  of  d  planes:  the  d-1  bounding  planes  of  A  and  the  plane  M.  Ltt  Rp  be  \ht 
plane  parallel  to  P  passing  through  Op.  We  note  4at  if  Rp  does  not  intersect  the  interior  of 
either  of  A  and  fl,  then  P  does  not  intersect  the  interior  of  at  least  one  of  these  regions.  We 
also  note  that  Rp  either  intersects  the  interior  of  both  A  and  fl,  or  it  intersects  the  interior  of 
neidier  A  nor  B.  In  addition,  these  two  observations  apply  in  any  affine  subspace  containing 
Op.  Let  C*-'  denote  the  subspace  M  H  ^^  H  '"  H  '"'"'"^  obtained  by  choosing  F^  to  be 
the  level  j  bounding  plane  of  A. 

Lemma  3:  For  each  i,  l^irsd-2,  there  is  a  set  Sj"',  such  that  the  following  holds.  For 
each  plane  P  in  S/"',  in  the  corresponding  subspace  Q''~',  the  corresponding  plane  Rp  does 
not  intersect  the  interior  of  either  of  the  regions  divided  by  P.  Furthermore,  Sf~_'^  C  5/^"'*^ 
(we  define  j^  =  0). 

Proof:  The  proof  is  by  induction  on  i.  For  the  base  case,  j=  1,  the  proof  is  similar  to  that  of 
lemma  1.  Let  /»  be  a  plane  in  S^.  The  subspace  C*"^  is  a  plane.  Consider  the  d- 1  bound- 
ing planes  for  A  and  fl.  With  respect  to  these  planes,  let  C  (resp.  D)  be  the  region  that  is  on 
the  opposite  side  as  A  (resp.  B)  of  the  level  7  plane,  for  l^j^d-1  and  j* d-1,  and  on  the 
same  side  for  ;  =  rf-l.  We  call  A,  fl,  C,  D,  the  regions  associated  with  P  in  Q'^'K  In  Q''~\ 
the  regions  A,  B,  C,  D,  are  the  four  regions  separated  by  the  following  two  lines:  the  inter- 
sections with  C*"^  of  the  level  d-2  and  level  d-1  bounding  planes.  The  origin  Op  is  the 
intersection  of  these  two  lines,  and  Rp  fl  Q'''^  »s  the  line  through  Op,  parallel  to  P  fi  Q''~^- 
We  observe  that  in  Q'''^,  Rp  does  not  intersect  the  interior  of  two  of  the  regions,  either  the 
mixed  or  the  unmixed  regions  (with  respect  to  the  level  d-2  and  the  level  d-1  planes).  Sup- 
pose that  it  does  not  intersect  the  mixed  (resp.  unmixed)  regions. 

Let  Q  be  any  plane  in  S^.  In  G'*"^  the  level  Jk  plane  bounding  the  four  regions  associ- 
ated with  Q  is  parallel  to  the  level  *  plane  bounding  Ae  four  regions  associated  with  P ,  for 
l^k^d-l.  Also,  in  Q''~^,  Rq  is  parallel  to  Rp.  Thus,  in  Q''~\  Rq  does  not  intersect  the 
interior  of  the  mixed  (resp.  unmixed)  regions  associated  with  Q.  We  choose  jj_i  accord- 
ingly. 

For  the  inductive  step,  suppose  that,  in  the  subspace  C*"'"^^,  for  eadi  plane  P  in  5^"'*^ 
the  corresponding  plane  Rp  does  not  intersect  the  interior  of  either  of  the  regions  divided  by 
/>.    Let  P  be  a  plane  in  Sj~''*'^.    Let  A  and  B  be  the  regions  divided  by  P.    The  region  C 

(resp.  D)  is  defined  to  be  the  region  with  the  following  signature:  the  j'-th  bit  is  opposite  to 
the  j-th  bit  of  A's  signature,  for  lsj<d-i,  and  is  the  same  for  d-i^j^d-1.   In  the  space 


Qi-i-*-i-  corresponding  to  P,  C  (resp.  D)  is  the  same  region  as  A  (resp.  B).  By  the  inductive 
hypothesis,  in  C''~'+\  Rp  does  not  intersect  the  interior  of  either  A  or  fl  (and  hence  of  either 
C  or  D).  We  claim  that  in  Q''~',  Rp  can  intersect  the  interior  of  the  regions  in  at  most  one  of 
the  pairs:  A  and  B,  C  and  D.  For  suppose  that,  in  Q''~',  Rp  intersects  the  interior  of  both  A 
and  C.  Consider  two  points  on  Rp  f)  C*"',  one  in  the  interior  of  A  and  one  in  the  interior 
of  C.  The  line  joining  these  two  points  intersects  g**"'"^^  at  a  point  in  the  interior  of  A 
(equaUy,  we  can  view  it  as  being  in  the  interior  of  C).  But  this  line  lies  on  Rp;  this  contrad- 
icts the  fact  Aat,  in  Q''~'*^,  Rp  does  not  intersect  the  interior  of  A.  Similarly,  in  C*"',  Rp 
intersects  the  interior  of  at  most  one  of  B  and  D.  Suppose  that,  in  Q*^'',  Rp  intersects  the 
interior  of  A.  we  argue  that  in  C**"',  Rp  also  intersects  the  interior  of  B.  For  consider  a 
point  in  C*"',  on  Rp,  in  the  interior  of  A ,  and  consider  the  line  joining  this  point  to  Op.  This 
line  passes  through  the  interior  of  B;  also  the  line  lies  on  Rp.  Hence,  in  G**"',  if  Rp  intersects 
the  interior  of  A,  then  in  C*"',  Rp  also  intersects  tfje  interior  of  B.  We  conclude  that,  in 
Q''~',  Rp  intersects  the  interior  of  regions  in  at  most  one  of  the  pairs  A  and  B,  C  and  D. 
Suppose  that  Rp  does  not  intersect  the  interior  of  the  mixed  (resp.  unmixed)  regions,  with 
respect  to  the  bounding  level  d-i-1  and  level  d-i  planes. 

Let  0  be  a  plane  in  Sf''""^.   In  C*"',  the  level  k  bounding  plane  for  the  four  regions 

associated  with  Q  is  parallel  to  the  level  k  bounding  plane  for  the  four  regions  associated  with 
/*,  for  l^k^d-1.  Also,  in  C*"',  Rq  is  paraUel  to  Rp.  Thus,  in  Q"'',  Rq  does  not  intersect 
the  interior  of  the  mixed  (resp.  unmixed)  regions.  We  choose  7j_,  accordingly,  a 

Corollary:  The  partition  is  avoiding. 

Proof:  Let  A/  be  a  plane.  We  show  that  the  partition  defines  at  most  2''-!  regions  on  M. 
We  can  then  deduce  that  M  intersects  at  most  2^-1  of  the  regions  defined  by  the  partition. 
Since  this  is  true  for  any  plane  M,  the  partition  is  separating. 

To  obtain  2''  regions  on  M,  we  would  need  for  each  level  d  plane  Q,  that  in  M,Q  inter- 
sects both  the  regions  divided  by  Q.   Consider  P,  the  plane  in  Sj .   By  lemma  3,  in  A/,  the 

corresponding  plane,  Rp,  does  not  intersect  either  of  the  regions  divides  by  P.  Hence,  in  M, 
P  intersects  at  most  one  of  the  two  regions  divided  by  i*.  The  result  follows.  □ 


3.  The  PatHdoD  b  Separating 

We  prove  that  by  a  judicious  choice  of  the  level  d  planes  we  can  obtain  an  n/1'^- 
partition.  Our  proof  will  proceed  by  giving  a  construction;  it  should  not  be  viewed  as  an 
algorithm,  however.  That  will  be  given  in  the  next  section.  There  are  two  major  concerns  in 
our  construction.  The  first  is  to  select  each  plane  in  such  a  way  that  it  is  unique.  The  second 
concern  arises  because  our  construction  will  involve  certain  parameters  that  may  vary.  The 
concern  is  to  ensure  that  the  planes  vary  continuously  with  these  parameters,  in  a  sense  that 
is  made  precise  below. 

We  build  the  partition  by  constructing,  in  turn,  the  level  1  planes,  the  level  2  planes,  ..., 
the  level  d  planes.  Since  constructing  the  level  /  planes,  for  i<d,  is  a  lower  dimensional 
problem,  it  suffices  to  describe  how  to  construct  the  level  d  planes.  We  will  show  that  given 
an  orthogonal  axis  system  x^,  ...,  x^,  a  J-dimensional  parallel  planes  partition  exists.  In  addi- 
tion, we  will  show  that  the  partition  is  unique,  in  the  following  sense.  We  discard  those 
points  lying  on  level  i  planes,  for  i<d.  Let  P  be  a  level  d  plane  being  constructed.  Suppose 
P  is  intended  to  bisect  point  sets  of  size  m^  and  mj,  contained  in  regions  A;  and  A^,  respec- 
tively. Then  we  require  that  P  pass  through  the  [(m,-i-l)/2j-th  point,  in  region  A,,  j  =  1,2 
(ordering  the  points  from  left  to  right  with  respect  to  F).  It  is  convenient  to  refer  to  a  plane 
P,  that  divides  the  point  sets  in  this  way,  as  bisecting  the  point  sets.  It  is  easy  to  see  that  this 
construction  provides  an  n/2''-partition. 

For  each  level  d  plane  P  we  will  define  a  series  of  angles  6^_i,  9j_2,  ...,  6;,  that 
describe  the  plane's  orientation.  (Strictly  speaking,  we  should  also  index  the  angles  8,  by  the 
index  for  the  corresponding  plane;  as  no  ambiguity  results,  we  omit  this.)  6^-|  and  6^_2  will 
be  Ae  same  for  all  planes.  For  each  »,  Is  j^  J-3,  and  for  each  j,  Osy  <2',  ^i-i-i  will  be 
the  same  for  all  planes  in  5y  ~'.  This  restriction  will  correspond  to  the  planes  in  5^''  being 
paraUel  in  the  subspace  /"^  H  ■"  H  P''~'~^-  We  define  the  angles  one  by  one. 

We  start  by  defining  9^_i.  Consider  the  plane  R''~^  =  ^^  D  "■  D  ^''~^-  (See  figure 
3).   In  R'''^,  let  y^.^  be  the  normal  to  P^~^;  define  8^_i  to  be  the  angle  between  Y^.^  and 

We        define        0^_,,        2^i^d—l,        similarly.  Consider        the        plane 

R"-'  =  P^  n  •••  n  '"'"'"^  n  ^"''^^  n  •■•  n  ^"^  where  we  choose  P^,  for  Is  its  J- 1, 
to  be  the  level  ik  plane  bounding  both  of  the  regions  divided  by  P.  (For  the  regions  divided 
by  planes  in  Sf~'*^,  all  the  level  it  bounding  planes  are  parallel  in  the  subspace 
P^  n  •••  n  P"'''^-)    In  H'''',  let  Yj.,  be  the  normal  to  P''-';  define  e^_,  to  be  the  angle 


between  y^.,  and  P^. 

We  need  to  show  that  the  angles  6,,  Is/s  J-1,  produce  the  parallelism  we  require. 

Lemma  4:  If  9^.^  and  e^_2  are  the  same  for  all  level  d  planes,  and  if,  for  each  »,  Is/sJt, 
for  eadi ;,  0s;<2',  e^_,_2  is  the  same  for  all  planes  in  Sj~\  then,  for  each  A,  OsA<2*,  the 
planes  in  5^~*  are  parallel  in  P^  f)  •••  fj  i»''"*~^  for  Osts<f-3. 

Proof:  We  prove  the  result  by  induction  on  it.  The  result  is  clearly  true  for  ik=0.  So  sup- 
pose is  is  true  for  ifc= «'- 1;  we  show  that  it  is  also  true  for  t=  i.  But  this  is  immediate  by  the 
construction  of  ^^-i-i.  n 

Let  6',  Isis  J-1,  denote  the  angles  producing  an  n/2''-partition.  We  wish  to  show 
that  the  angles  SJ-i  and  62-2  ^^i^^  ^^^  "^  unique,  and  likewise  for  each  angle  6j_/_2'  associ- 
ated with  a  set  S/"',  for  0sj<2'  and  lsjsd-3.  Actually,  we  will  show  that  each  term 
tand*  (denoted  r*)  exists  and  is  unique;  but  this  is  the  same  result.  We  extend  the  notation  to 
make  t'  a  function  of  tan8y  (denoted  fy),  for  }>i,  as  follows.  For  each  set  5^*\  given  the 
corresponding  values  Cy,  for  j>i,  t'  denotes  the  value  of  t,  producing  the  required  parallelism 
for  the  planes  in  S'^'*'^.  Before  proving  the  existence  and  uniqueness  result  we  need  some 
definitions  and  a  technical  lenama. 

Definition:  Let/(x;t,...,xJ  be  a  fimction  of  k  variables.  /  is  unbounded  in  x,  if,  when  x,  tends 
to  »,  and  Xj  is  bounded,  for  j^i,  then /  tends  to  one  of  ±<»;  and  when  X[  tends  to  -«,  /  tends 
to  one  of  ±ao.  /  is  a  strict  function,  if,  as  a  function  of  x,,  /  is  either  strictly  increasing  or 
strictly  decreasing,  for  l^t'^Jt.  (f  could  be  strictly  increasing  in  one  variable  and  strictly 
decreasing  in  another  variable.)  If  /  is  a  strict  function,  the  signature  of/  is  defined  to  be  the 
ik-bit  number,  whose  t-th  bit  (from  the  right)  is  a  one  if  and  only  if  /  is  an  increasing  function 
of  x„  for  Is  is*. 

Lemma  5:   Let  the  function  f(xi^ xJ  be  strict  and  have  signature  a.    Also,  suppose/  is 

unbounded  and  continuous  in  each  variable.  Consider  the  fimction  gix^,  .  .  .  ,X2),  which, 
when  defined,  satisfies  /(x4,...,X2,^)  =  0.  Then  g  is  defined  for  all  x^,  ...rX2'  •*  unbounded 
and  continuous  in  each  variable,  and  is  a  strict  function.  If  a= 2t,  f  has  signature  t;  while  if 
a  =  2r+l,  g  has  signature  t,  the  complement  of  t. 

Proof:  We  show  that  ;  is  a  continuous  function  of  x/,  for  2^i^k;  the  other  claims  are 
straightforward  and  are  left  to  the  reader.  Without  loss  of  generality  suppose  that  /  is  an 
increasing  fimction  of  X|  and  x^  (hence  ;  is  a  decreasing  function  of  x,).  We  show  that  j  is  a 
continuous  function  of  x,  from  below;  the  proof  of  continuity  from  above  is  identical.  Let 
g~  =  lim  g(x^,...,x,~h,...,X2).   And  let  g  =  g(x^ X2).  Since  ^  is  a  decreasing  function  of 
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'/.  8~^8-  By  the  continuity  of/,  ]imf(x^,...^i-h,...,X2,g~)  =  fix^,...,x,,...,X2,g~).  Since 
/is  an  increasing  function  of  zi,/(ii,...,jc„. ..,X2,j~)  i/(Xi,...,ac„...,JC2,«)  =  0.  Since  g  is  a 
strictly  decreasing  function  of  x„  and  /  is  a  strictly  increasing  function  of  x^, 
fix^,...,x,-h,...,X2,g-)<fixt,...,x,-t,...,X2,gixt,...,x,-h,...,X2))  =  0.  Letting  8^,  we 
obtain,  by  the  continuity  of  /,  f(xi„...,x„...,X2,g~)  ^  0.  We  deduce  that 
fixt,...,x„...,X2,g-)  =  0,  and  hence  g-  =  g.  d 

For  each  level  d  plane  P,  and  the  associated  angles  9*,  ...,  e^-i,  we  will  show,  by 
repeated  use  of  lemma  5,  that  t'  is  a  continuous,  strict,  unbounded  function  of  tj,  for  j>i. 

We  start  by  considering,  for  each  region  t/^"^  the  projected  point  set,  obtained  by  pro- 
jecting the  point  set  contained  in  if^~^,  in  turn,  in  the  direction  e**"^  onto  P''"^  ,..,  in  the 
direction  e**"'  onto  i"*"',  ...,  in  the  direction  6^  onto  P^.  Let  x„^a  b*  ^^  median  of  the  pro- 
jected point  set,  defined  as  follows:  let  the  set  of  projected  points  comprise  m  points;  then 
x„^    is    required    to    be    the    [(m  +  l)/2j-th    point   from    the    left.     We   will    show    that 

^mtdi'd-i O  '*  continuous,  strict,  and  unbounded,  as  a  function  of  each  f,,  Isisd-l 

(where  we  consider  the  variables  r,  associated  with  the  plane  P  bisecting  U''~^). 

We  define  (with  respect  to  P)  the  origin  to  be  some  point  in  the  intersection  of  the 
planes  bounding  the  region  U''~^  (this  intersection  is  a  line  parallel  to  P^  f)  ■■■  (~)  P''~^).   It 
is  useful  to  define  a  new  set  of  (non-orthogonal)  coordinate  axes  (again  with  respect  to  P), 
namely,  y^,  ...,  y^.  yt  is  parallel  to  the  line  Y'  defined  above,  for  l:sisd-l.  y^  is  parallel  to 
^^  n    ■■  n  f"''^-    All  the  axes  pass  through  the  origin,  and  take  the  origin  as  their  zero. 
Also,  for  Is  jsi-1,  axis  >,  is  oriented  to  match  the  orientation  of  the  level  i  plane,  bound- 
ing the  region  V''^  at  hand,  and  yj  is  oriented  to  match  the  orientation  of  P.   Consider  a 
point  (ai,...,aj)  in  U''~^,  the  coordinates  being  with  respect  to  the  axes  y,,  Isjsrf.   (These 
coordinates  are  independent  of  the  angles  8^,  ...,  S^.^.)   Under  the  projection  described  in 
the  last  paragraph  the  point  is  mapped  to  the  projected  point  aji+  ■  ■  •  +fld-i^<<-i+fld-   We 
note  that  each  point  in  U''~^  has  the  same  signature.  Thus  increasing  some  r,  has  the  effect  of 
increasing  or  decreasing  every  projected  point  from  U^~^;  we  conclude  that  x„^(tfi.^,...,tj  is 
a  strict  function  and  that  it  has  the  same  signature  as  the  set  U^~^.   Also,  it  is  dear  that  x„^ 
is  a  continuous  function  of  each  tj,  l^J^d-1.  Suppose  that  l/''^  lies  on  the  positive  (resp. 
negative)  side  of  its  bounding  level  j  plane.   If  we  bound  every  variable  except  for  tj,  and  let 
this  latter  variable  tend  to  »,  then  x„^  tends  to  <*>  (resp.  -•).   Similarly,  if  we  let  tj  tend  to  -» 
then  x„^  tends  to  -«  (resp.  <*>).    Thus  x„^,  as  a  function  of  tj,  is  continuous,  strict,  and 
unbounded,  for  Irsj^d-l,  and  it  has  the  same  signature  as  the  region  U''~^. 
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We  can  now  obtain  the  following  result. 

Lemma  6:  Let  i<d-l,  r^.^ r,+2>  '^^  *  ^^^  ^r*^  ^  given.   Let  a{r)  denote  the  signature 

ofS'/K  Then: 

(1)  Si^jpose  that  r,+i  is  also  given.  There  is  one  choice  of  t„  denoted  r/  (resp.  f,^),  that 
achieves  the  required  parallelism  and  bisection  for  the  planes  in  5^^^  (resp.  52r+i)-  As 
functions  of  t,,  for  j>i,  t}  and  if  are  continuous,  strict,  and  unbounded.  In  addition,  r/ 
has  signature  (r(2r),  while  t}  has  signature  a(2r+l),  the  complement  of  a(2r). 

(2)  There  is  a  unique  choice  of  r,+i  (denoted  r/Vj  for  which  t}  =  tf,  and  further,  as  a 
function  of  tj,  j>i+l,  t'+i  is  continuous,  strict,  and  unbounded.  In  addition,  t'^^  has 
signature  <T{r). 

Proof:  We  prove  the  result  by  induction  on  j. 

The  base  case  will  be  for  (  =  1.  We  prove  each  of  the  claims  in  turn. 

Claim  1:  Let  P  be  the  plane  in  5^,.  Let  A^  and  A^  be  the  two  regions  bisected  by  P,  with  A^ 
on  the  negative  side  of  its  bounding  level  1  plane.  Recall  that  the  signatures  of  A^  and  A^  are 
complementary.  Let  x^,^  be  the  value  of  x„^  (as  defined  above)  corresponding  to  region  A,, 
1  =  1,2.  As  shown  above,  as  a  function  of  tj,  for  lsjsd-1,  x,'„,^  is  continuous,  strict,  and 
unbounded;  also  it  has  the  same  signature  as  region  A,.  Hence  x^,^  and  x^,^  have  comple- 
mentary signatures.  Define  diff^  =  xj;^-xl^;  as  a  function  of  tj,  for  l^j^d-l,  diff^  is 
continuous,  strict,  and  unbounded;  also,  diff^  has  the  same  signature  as  jcj,^. 

When  diff^  =  0,  P  bisects  the  point  sets  in  each  of  the  regions  A,,  i  =  1,2.  rf  is  the  value 
of  f^  for  which  diff^  =  0;  by  lemma  5,  t{  is  a  continuous,  strict,  and  unbounded  function  of  tj, 
for  7>1.  We  recall  that,  by  definition,  the  signature  a(2r)  of  S],  consists  of  the  first  d-2 
bits  (from  the  left)  of  the  signature  of  A^.  (See  the  remark,  following  lemma  2,  defining  sig- 
natures  for  sets  5'*^^.)  In  addition,  we  note  that  the  rightmost  bit  of  the  signature  of  A^  is  a 
zero  (since  it  is  on  the  negative  side  of  its  level  1  bounding  plane);  thus  A^  has  signature 
a(2r)0  (=2a(2r)).  As  observed  above,  just  before  the  current  lemma,  A^  and  x^  have  the 
same  signature,  as  does  diffQ.  So,  by  lemma  5,  t{  has  signature  a(2r). 

We  conclude  that  i{  is  a  continuous,  strict,  unbounded  function  of  tj,  for  ]>!,  with  sig- 
nature a(2r).  Similarly,  we  can  prove  that  r^  is  a  continuous,  strict,  unbounded  function  of 
tj,  for  j>l,  with  signature  (T{2r+1). 

We  show  that  a(2r-l-l)  is  the  complement  of  a(2r).  Let  A^  have  signature  tOO  (recall 
that  the  plane  in  5|^  divides  unmixed  regions  with  respect  to  the  level  2  and  level  1  bounding 
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planes,  and  that  we  assumed  A^  was  on  the  negative  side  of  its  bounding  level  1  plane).  So 
a(2r)  =  tO.  ^2  has  signature  fll  and  the  regions  divided  by  the  plane  in  Sl^+^  have  signa- 
tures tOI  and  tIO.  So  a(2r+l)  =  t1;  this  is  the  complement  of  tO  =  <T(2r),  as  claimed. 
Claim  2:  Consider  diff^  =  ^i-'i-  As  a  function  of  f,,  for  j>l,  diff^  is  continuous,  strict  and 
unbounded.  Also,  diff^  has  signature  <y(2r).  Hence,  by  lemma  5,  rj,  the  value  for  which 
diff^i  =  0,  as  a  fimction  of  tj,  for  J>2,  is  continuous,  strict,  and  unboimded.  Since 
a(2r)  =  tO,  we  deduce,  by  lemma  5,  that  fj  lias  signature  t.  By  definition,  the  signature 
a{r)  of  S^  consists  of  the  first  d-3  bits  of  the  signature  of  A^;  thus  a(r)  =  t. 

For  the  inductive  step,  suppose  that  the  result  is  true  for  j  =  *-l.   We  show  that  it  is 
also  true  for  i  =  it.  We  prove  the  two  claims  in  turn. 

Claim  1 :  Almost  all  of  this  claim  has  been  proved  in  claim  2  of  the  inductive  hypothesis  for 
Jk- 1.  It  only  remains  to  prove  that  a(2r+ 1)  is  the  complement  of  a(2r).  Among  the  regions 
divided  by  planes  in  Sit^  let  A  be  the  region  on  the  negative  side  of  its  level  )  bounding 
planes,  for  ls;SJk.  The  first  d-k-\  bits  (from  the  left)  of  the  signature  of  A  comprise  Ae 
signature  of  S^*^  We  note  that  A  has  the  (J-l)-bit  signature  tO**\  and  hence  a(2r)  =  T0. 
The  regions  divided  by  planes  in  S^t^i  indude  the  region  with  signature  tIO*.  Thus 
a(2r+l)  =  Tl,  the  complement  of  a(2r). 

Claim  2:  To  show  that  rj^+i  is  a  continuous,  strict  and  unbounded  function,  we  argue  as 
above,  replacing  diffi  by  diffi,  =  tj;-tl.  We  can  also  conclude  that  t^+i  has  signature  t;  it 
remains  to  show  that  a{r)  =  t.  As  argued  in  claim  1,  region  A  has  signature  tO*'^^  Since  A 
is  on  the  negative  side  of  its  bounding  level  j  planes,  for  Isjsjt,  a(2r)  =  T0.  Also,  A  is  on 
the  negative  side  of  its  boimding  level  j  planes  for  Isj^k+l,  and  A  is  a  region  bisected  by 
a  plane  in  S*"^^.  Thus  a(r)  =  r.  d 

Lemma  6  shows  that  for  each  level  d  plane,  the  corresponding  values  f',  fi,  ...,  fj-i 
exist  and  are  unique.  We  deduce  that  each  level  d  plane  exists  and  is  unique,  and  hence  a 
separating  parallel  planes  partition  exists  and  is  unique,  given  the  coordinate  axes. 

4.  Algorithms  for  partitioning 

We  start  by  describing  a  deterministic  algorithm  in  2  dimensions.  It  will  be  used  to 
motivate  and  illustrate  the  algorithm  in  higher  dimensions.  This  is  not  our  most  efficient 
algorithm.  Later  we  will  show  how  to  achieve  greater  efficiency. 

Recall  the  problem.  We  have  two  sets,  each  comprising  at  most  n/2  points,  separated 
by  the  line  x  =  0.   We  seek  a  line  L  that  simultaneously  bisects  both  these  sets.   Suppose  this 
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line  makes  an  angle  6*  with  the  line  y  =  0.    We  show  how  to  find  6*.    Or  rather,  we  show 
bow  to  find  tane*,  denoted  t'.  It  is  then  easy  to  detennine  L. 

Our  method  is  to  project  the  points,  parallel  to  L  (currently  unknown),  onto  the  line 
x  =  0.  We  then  determine  the  median  of  each  set  of  projected  points  (by  sorting  the  points); 
this  computation,  as  a  side  effect,  will  yield  9*. 

Recall  that  the  point  (a,b),  under  projection,  parallel  to  L,  is  mapped  to  at'  +  b.  The 
problem  is  to  sort  expressions  of  the  form  at'  +  b,  when  we  do  not  know  the  value  of  /*. 
This  problem  can  be  solved  using  the  methodology  developed  by  Megiddo  for  solving  such 
parameterized  problems  [Ml].  Rather  than  explain  his  general  methodology,  we  instead  con- 
centrate on  a  restricted  case  to  which  a  methodology  discovered  by  Cole  [C2]  applies.  The 
latter  applies  to  sorting  problems  satisfying: 

(1)  Comparisons  can  be  ordered.  That  is,  given  a  series  of  comparisons  C^,  ...,  C^, 
there  exists  an  ordering  of  these  comparisons,  denoted  C^(i,sC^(2)^---sC„(p),  where 
the  ordering  has  the  following  meaning.  We  let  Dj  denote  C„(,^,  for  l^i^p.  Suppose 
D,  is  the  comparison  yn'.ya-  K  the  result  of  comparison  D,  is  that  >/i<>q  (resp.  yii^yi2) 
then  the  resxilt  of  comparisons  Cj,  for  j<i,  is  that  yji<yj2  (resp.  yji^yj2)',  conversely,  if 
yti>yt2  (resp.  yn^ya),  then  Ae  result  of  comparisons  Cj,  for  j>i,  is  Aatyj^>yj2  (resp. 

(2)  The  median  of  C^,  ...,  C„,  in  the  ordering  described  in  (1),  can  be  computed  in 
M{n)  time. 

(3)  The  result  of  a  comparison  can  be  computed  in  C(n)  time. 

For  sudi  sorting  problems,  we  can  carry  out  the  sort  in  0(nlogn  +  \ogn[M(n)  +  C(n)]) 
time.  (In  [C2]  it  is  implicitly  assumed  that  M(n)  =  0(n);  this  will  not  always  be  the  case 
here.)  We  refer  to  this  algorithm  as  the  parameterized  sorting  algorithm. 

We  apply  this  sorting  algorithm  to  o\ir  problem.  That  is,  we  show  that  the  problem  can 
be  put  in  a  form  satisfying  (1),  and  we  give  algorithms  for  (2)  and  (3). 

First,  we  explain  what  it  means  to  compare  two  projected  points  Pi  =  aj'-^b^  and 
Pi  ~  02t''^b2.  This  is  simply  a  matter  of  deciding  which  of  p^  and  p2  is  Ae  larger;  that  is, 
we  need  to  determine  the  sign  of  ct'+d,  where  c  =  a^-  02  a°d  d=b^-  62-  If  c  =  0  this  is  sim- 
ply the  sign  of  d;  otherwise,  we  might  as  well  compute  the  sign  of  f*  +  *,  where  e  =  d/c. 

We  order  those  comparisons,  not  yet  resolved,  according  to  the  ordering  of  the 
corresponding  terms  *,  considered  as  numeric  values.  Next,  we  describe  how  to  find  the 
median  of  n  comparisons,  that  is  of  n  expressions  of  the  form  t'  +  e,,  Isfsn.   This  can  be 
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done  in  linear  time  by  computing  the  median  of  the  numbers  «;. 

The  final  step  is  to  evaluate  the  sign  of  an  expression  t'  +  e.  We  choose  t  to  satisfy 
J^g  =  0  (let  e  =  7).  We  project  each  set  of  points  under  the  angle  0,  and  compute  the  median 
of  each  set  of  projected  points,  in  linear  time.  According  to  the  relative  ordering  of  the  two 
medians,  we  determine  that  6'  is  larger  or  smaller  than  9  (see  lemma  6),  and  hence  we  deter- 
mine the  sign  of  t'  +  e.  (If  the  two  medians  are  equal,  we  have  found  the  partitioning  line  L, 
and  we  stop  the  algorithm  at  this  point.)  Notice  that  this  computation  has  reduced  the  range 
in  which  6*  can  lie  (by  determining  one  of:  e*<e,  e*>e).  We  keep  track  of  this  range; 
when  the  sort  is  complete,  the  range  for  0*  will  have  been  reduced  to  a  single  value  (since, 
when  the  sort  is  completed,  the  medians  for  the  two  sets  of  points  are  shown  to  be  equal). 

We  have  shown  that  we  can  sort  the  projected  points  using  the  parameterized  sorting 
algorithm.  We  deduce  that  we  can  sort  the  projected  points  in  0(nlog«)  time.  We  have  also 
foimd  the  partitioning  line:  it  is  the  line  passing  through  the  medians  of  the  two  sets  of 
points. 

Remark:  There  is  a  linear  time  algorithm  for  this  problem  due  to  Megiddo  [M2].  We  will 
incorporate  this  algorithm  into  our  solution  later;  the  above  algorithm  is  intended  to  help  the 
reader  understand  the  general  algorithm  given  below. 

We  turn  to  the  general  problem  in  d  dimensions.  We  recall  there  are  2''"^  regions 
t/^~^.  With  each  such  region  we  associate  the  sequence  of  angles  6',  ej,  ...,  SJ-i,  associated 
wiA  the  level  d  plane  dividing  the  region.  For  each  region,  we  project  the  points  in  the 
region  under  the  corresponding  d—  1  angles,  obtaining  a  set  of  projected  points.  For  each  of 
these  l**"^  sets  we  want  to  sort  the  projected  points,  and  thereby  determine  their  medians  (as 
in  the  2-dimensional  case,  the  sort  will  also  yield  the  actual  values  of  6*,  im^d-1).  We 
recall  that  the  point  sets  contained  in  regions  with  opposite  signature  should  have  equal  medi- 
ans, that  for  point  sets  bisected  by  planes  in  5]"*"^  the  angle  6'  is  the  same,  for  l^i^d-2  and 
0^j<2'^~'~^,  and  that  SJ-i  is  the  same  for  all  the  point  sets.  These  constraints  will  be  used 
to  guide  the  sorting  algorithm. 

Under  projection,  the  point  (ai,...,aj)  is  mapped  to  ajl+  •  •  •  fld-i'j-i  +  flj-  Thus 
comparing  two  projected  points  is  a  matter  of  determining  the  sign  of  an  expression  of  the 
form  bjl+  ■  ■  •  b^^j2-i'^^d-  ^^ '  ^  *^  smallest  index  such  that  b,  is  non-zero;  equally, 
we  could  compare  the  two  projected  points  by  determining  the  sign  of 
r;4-c,+ir,%i  +  ...  +  Crf_ir;_i-t-c^,  where  Cy  =  b/6,. 

At  this  point,  rather  than  continue  explaining  how  to  sort  the  projected  points,  we 
describe  two  more  general  problems,  and  show  how  to  solve  them  by  mutual  recursion.   We 
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will  show  that  the  sorting  problem  is  an  instance  of  the  general  problem. 

The  first  problem,  denoted  S(j',ii,...,i^),  is  the  problem  of  finding  the  median  of  (usu- 
ally by  sorting)  the  n  expressions  ay«';+  •  •  •  +a/^-i'rf-i+flw.  l^i^n,  under  the  following 

d-i 

constraints.   For  i*i^,  l^r^k,  t(  =  t',  while  for  1  =  1,.,  r,  =  'Z^u',^  IsrsJt,  where  the  con- 

s>l 

straints  are  applied  recursively.  (Recall  that  for  eadi  set  5^"^^  8 J  is  defined  to  be  the  angle 
producing  the  required  parallelism  for  the  planes  in  S^"*"-^,  given  the  values  t,,  for  s>j.)  We 
require  that  i{>  ■  •  ■  >i^,  i^^j,  and  l^j^d. 

The  second  problem,  denoted  i?0;'i. •••.»»).  »s  the  problem  of  computing  the  sign  of  the 
expression  fly<}+  •  •  •  +0^-1^^-1+0^,  under  the  same  constraints  as  for  the  first  problem, 
and  the  same  requirements  except  that  here  iji>j. 

We  solve  each  problem  in  turn.  For  the  first  problem,  assuming  ii^>j  and  j<d,  we 
apply  the  parameterized  sorting  algorithm.  (We  discuss  the  cases  1^  =  ]  andj  =  d  below.)  A 
comparison  of  two  expressions  consists  of  determining  the  sign  of  an  expression  of  the  form 
bjtj-¥  •  •  ■  +6^-i^rf_i+fcd-  Let  Jt  be  the  least  index  such  that  b^  is  non-zero.  Determining 
this  sign  is  equivalent  to  computing  the  sign  of  ^4  +  ^4+1/4+1+  •  •  •  +c^-i^<i-i+c^,  where 
Ci^bf/b^.  Since  resolving  a  comparison  amounts  to  finding  out  the  sign  of  this  expression, 
we  can  order  the  comparisons  by  the  numeric  order  of  the  corresponding  expressions. 
According  to  the  methodology  given  above,  we  have  two  tasks  to  carry  out.  The  first  is  to 
find  the  median  of  n  such  comparisons;  the  second  is  to  determine  the  sign  of  one  of  these 
comparisons  (the  median  one,  in  fact).  Having  carried  out  these  two  tasks  we  are  able  to 
resolve  one  half  of  the  comparison.  We  will  proceed  differently;  at  the  end  of  the  process  we 
will  still  be  able  to  resolve  one  half  of  the  comparisons.  In  fact,  all  we  need  for  the 
parameterized  algorithm  is  that  we  resolve  one  half  of  the  comparisons;  the  actual  method 
used  is  of  no  consequence  (see  [C2]).  The  method  we  use  is  to  divide  the  comparisons  into 
d—j  sets  (for  reasons  that  will  become  clear)  and  determine  the  sign  of  the  median  com- 
parison in  each  set;  we  can  then  resolve  one  half  of  the  comparisons  in  each  set,  and  hence 
we  resolve  one  half  of  the  comparisons  overaU. 

To  find  the  median  of  n  such  comparisons  we  would  sort  the  corresponding  expressions. 
If  the  leading  term  in  each  expression  were  tj,  we  would  solve  this  problem  by  sorting  the 
remainder  of  each  expression,  that  is  expressions  of  the  form  Cy+iry+^-i-  •  •  •  +c^-i'^-i+c^i 
under  the  same  constraints  as  before;  this  is  just  the  problem  50  + l;»i. •••>'*)•  But  the  lead- 
ing term  in  some  expressions  need  not  be  r^.  So  we  divide  the  expressions,  according  to  their 
leading  term,  into  d-j  sets.   For  each  of  these  sets  we  perform  a  sort,  and  then  we  evaluate 
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the  sign  of  the  median  expression  in  that  set.  (For  the  purposes  of  analyzing  the  ninning 
time  we  will  consider  each  of  these  sorts  to  be  an  instance  of  the  problem  5(/  +  l;»i.--.«i). 
and  each  of  the  sign  computations  to  be  an  instance  of  J?0;'ii-*<it)) 

U  j  =  d,  5((i;)  is  the  problem  of  fmding  the  median  of  n  values  a^^,  l^i^d.  We  can 
compute  this  in  linear  time.   If  ii=y,  we  substitute  for  r,  and  SO' ;'i> •••.'*)  becomes  the  prob- 

lMn50  +  l;ti,...,»i-i)- 

For  the  second  problem,  we  seek  the  sign  of  an  expression  ajtj+  •  •  •  +a^_ir^_i+a^, 
under  the  constraints  given  above.  Let  tj  be  the  value  of  tj  for  which  the  above  expression  is 
zero.  Suppose  tj  corresponds  to  the  set  Sl*^.  We  seek  to  find  the  values  tj^^  and  tj_^ 
corresponding  to  the  sets  S{r  and  S^^+a  when  tj  =  Tj.  (see  lemma  6).  If  tj^i  =  tj.-^  we  know 
fy=ry.  Otherwise,  if  tj.^-ij.i>Q  (resp.  <0),  we  know  tj-tj>0  (resp.  <0).  This  then 
determines  the  sign  of  the  expression.  To  find  tj^i  and  tj_^  we  compute  the  median  of  each 
of  2J  sets  of  projected  points,  each  such  set  comprising  the  points  in  one  of  the  regions 
bisected  by  a  plane  in  Sl"*"^,  under  the  constraints  already  mentioned,  except  that  we  replace 
the  constraint  tj  =  tj  by  the  constraint  ajtj+  ■  ■  •  -^0^-^^-1+0^  =  0.  But  this  is  2J  instances 
of  the  problem  S{l;ii,...,i^,j).  On  completing  all  the  sorts  we  will  have  obtained  the  values 
6j_|  and  8^_;,  and  thus  we  can  deduce  the  sign  of  the  expression  ajtj+  ■  •  •  -^Od-Jd-i'^'^d' 
under  the  constraints  given  above. 

The  original  sorting  problem  we  mentioned  can  now  be  seen  to  be  2''"^  instances  of 
5(1;).  It  remains  to  show  that  the  recursion  for  S  and  R  is  finite,  and  to  analyze  the  r\mrung 
time.  We  carry  these  out  together. 

Consider  50*;«ii-".«4);  if  'i-*J  andj<d,  it  induces  0([J-j]logn)  recursive  calls  to  each 
of  50'  +  l;»i..-..O  *^*^  ^0;'i>-"»*i)'»  while  if  i^=j  it  induces  one  recursive  call  to 
SO"+l;«i."-i»i+i)-  Finally,  if  j  =  d,  the  procedure  nms  in  linear  time.  And  consider 
l?(;;»i,...,i\);  it  induces  2^  recursive  calls  to  S{l;ii,...,ii^,j).  Consider  the  recursion  tree, 
where  we  condense  every  node  of  degree  one  and  every  node  labeled  R  into  its  parent.  We 
obtain  a  tree  whose  nodes  have  degree  0(2''logn).  If  we  can  bound  the  height  of  the  tree  by 
h  we  can  deduce  the  algorithm  runs  in  time  0(n(clogn)''),  where  c  is  an  appropriate  con- 
stant, because  each  leaf  represents  linear  work.  Thus,  if  we  denote  by  P50";»i.---,»i)  the 
maximum  path  length  at  a  node  corresponding  to  the  recursive  call  S(j;ii,...,i^),  we  have 

for;<ii,     PS(j;i^ ij  =  1  +  max[P50  +  l;«i i*),  PS(l;«i ikJ)] 

and  if  j  =  ;i,     PSijUi '*)  =  ^SO  +  l;»i.  •  •  •  .«i-i) 

and  P5(J;)  =  0. 
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We  claim  that  the  length  of  the  longest  path  from  a  node  representing  •fO;'i>->^)  is 
2<'-i  -  2)  2 '     -  2'"^.  It  is  easy  to  check  the  correctness  of  this  formula  by  induction.  Thus 

r-l 

the  height  of  the  tree  for  5(1;)  is  g(d)  =  2'''^  -  1,  which  gives  us  a  running  time  of 
0(n(clogn)*^'')),  where  c  is  a  suitable  constant,  (c  =  bl'',  where  b  is  the  constant  that  arises 
in  the  parameterized  algorithm;  this  constant  b  is  independent  of  the  problem  at  hand). 

As  commented  above,  Megiddo  gives  a  linear  time  algorithm  for  the  2-dimensiona] 
problem  [M2].  We  can  use  this  algorithm  for  the  problem  S(l;<f-l,</-2,...,2),  thereby 
reducing  all  the  nmning  times  by  a  factor  of  0(}ogn).  Below,  we  show  that  we  can  trim  off 
more  factors  of  O(logn);  this  latter  improvement  is  more  complicated,  however. 

Before  obtaining  this  improvement,  we  remind  the  reader  of  Megiddo's  solution  to  Ae 
2-dimensional  problem.  To  keep  its  application  to  the  ti-dimensional  problem  clear,  rather 
than  solve  it  in  a  dual  space  as  Megiddo  does,  we  solve  it  in  the  original  space. 
Problem:  Given  n^  red  points  above  the  line  >'  =  0,  and  rij  black  points  below  the  line  y  =  0, 
find  the  bisecting  line  with  ^n^  red  points  and  Sn^  black  points  to  its  left,  and  with 
sn^-n^-l  red  points  and  Sn2-n4— 1  black  points  to  its  right  (n^^n^,  11^^112). 
Soladon:  We  give  a  linear  time  reducing  algorithm  that  transforms  the  input  problem  of  size 
fi  into  an  instance  of  the  same  problem  of  size  at  most  7n/8.  Hence  the  whole  problem  can 
be  solved  in  linear  time.  We  describe  the  reducing  algorithm  as  a  sequence  of  seven  steps. 

(1)  In  linear  time,  compute  the  median,  y„^,  of  the  3* -coordinates  of  the  points. 

(2)  Pair  the  points,  except  for  those  points  on  the  line  y=y,„^,  each  pair  consisting  of  one 
point  below  y  =  y^ts  ^^  °^^  po">*  above  >=>'„^^.  Leftover  points  not  on  y  =  y^^,  if  any,  are 
paired  with  points  on  y=y„^.  The  remaining  points  on  y=y„^,  if  any,  are  termed  the 
unpaired  points. 

(3)  Join  each  pair  of  points  by  a  line  and  compute  its  slope;  in  linear  time,  compute  s,  the 
median  slope. 

(4)  Determine  whether  the  bisecting  line  has  slope  greater  or  less  than  j,  as  follows.  In 
linear  time,  compute  the  line  L,  of  slope  s,  that  passes  through  the  (n3+l)-th  red  point  (the 
points  being  ordered  from  left  to  right  with  respect  to  L).  If  there  are  more  than  rti,  (resp. 
iij-n^-l)  black  points  to  the  left  (resp.  right)  of  L  then  the  dividing  line  has  slope  less 
(resp.  greater)  than  s.  (See  figure  4.)  The  only  other  possibility  is  that  L  is  the  bisecting  line 
sought,  in  which  case  the  algorithm  terminates. 


19 


(5)  Without  loss  of  generality,  suppose  that  the  bisecting  line  has  slope  greater  than  s.  Con- 
sider the  lines,  constructed  in  (3),  of  slope  less  than  *.  For  each  such  line  compute  the  point 
at  which  it  crosses  the  \iney=y„^.  To  these  points  add  the  unpaired  points  from  (2),  form- 
ing a  set  5  of  points.  In  linear  time,  compute  the  median,  x„^,  of  the  x-coordinates  of  the 
oointsinS.  Letp  =  (x^^,>„^). 

(6)  Determine  if  the  bisecting  line  passes  to  the  left  or  the  right  of  p ,  as  follows.  Suppose 
that>^,stO  (resp.  sO).  Determine  the  line  K,  through />,  also  passing  through  the  (n^-t-l)- 
th  black  (resp.  (n3-t-l)-th  red)  point,  as  follows.  For  each  black  (resp.  red)  point,  compute 
the  slope  of  the  line  joining  this  point  to  p.  Compute  the  (n4-i-l)-th  (resp.  (n3-t-l)-th)  slope, 
j'.  IT  is  the  line  of  slope  s'  passing  through  p.  For  each  red  (resp.  black)  point,  determine 
to  which  side  of  K  it  lies.  If  more  than  n^  red  (resp.  n^  black)  points  lie  to  the  left  of  K  then 
the  bisecting  line  is  to  the  left  of  p,  while  if  more  than  n^-n^- 1  red  (resp.  n2-n4-  1  black) 
point;^  lie  to  the  right  of  K  then  the  bisecting  line  is  to  the  right  of  p.  (See  figure  5.)  The 
only  other  possibility  is  that  K  is  the  bisecting  line  sought,  in  which  case  the  algorithm  ter- 
minates. 

(7)  Without  loss  of  generality,  suppose  that  the  bisecting  line  passes  to  the  left  of  p.  Con- 
sider the  lines  constructed  in  (3),  of  slope  less  than  5,  to  the  right  of  p.  Each  such  line  is 
defined  by  two  points.  The  endpoint  on  or  above  y=y,^  lies  to  the  right  of  the  bisecting 
line.  (See  figxire  6.)  The  same  is  true  of  those  points  on  y=y„^  to  the  right  of  p .  There  are 
at  least  n/8  such  points.  Since  we  know  to  which  side  of  the  bisecting  line  these  points  lie  we 
can  eliminate  them,  thereby  producing  a  smaller  instance  of  the  same  problem.  Similarly,  in 
the  other  cases,  we  also  can  eliminate  at  least  n/8  points,  n^,  112,  n^,  and  n^  are  adjusted 
accordingly. 

We  describe  the  solution  in  3  dimensions  next;  it  is  similar.  We  apply  the  above  algo- 
rithm twice:  once  to  the  points  in  the  two  regions  divided  by  the  plane  in  Sq  and  once  to  the 
points  in  the  two  regions  divided  by  the  plane  in  Sf.  We  project  from  3  dimensions  to  2 
dimensions  by  mapping  the  points  parallel  to  62:  a  point  {0^,02, a^)  is  mapped  to  (0^,0212"^ a^) 
(which  we  treat  as  the  y  and  x  coordinates,  respectively).  We  carry  out  steps  1  through  7,  as 
given  above,  but  with  the  points  being  parameterized  by  r^  here;  this  necessitates  some 
changes,  described  below. 

We  face  two  difficulties  in  using  the  above  algorithm.  One  is  to  compute  the  k-\h  larg- 
est of  terms  involving  rj.  We  solve  this  by  using  the  algorithm,  given  above,  for  5(2;)  in  3 
dimensions;  this  algorithm  takes  time  O(nlogn).  A  second  difficulty  is  to  evaluate  the  sign  of 
an  expression  t^  +  b.    We  solve  this  by  using  the  algorithm,  given  above,  for  /f(2;)  in  3- 


dimensions;  this  algorithm  takes  time  0(n).  We  refer  to  these  two  procedures  as  computa- 
tions  1  and  2,  respectively. 

Next,  we  describe  the  changes  that  need  to  be  made  to  each  of  the  seven  steps  of 
Megiddo's  algorithm. 

Steps  1  and  2  are  the  same  as  in  the  2-dimensional  algorithm. 

Step  3.  Let  {a^,aQ  +  a2tl),  (^i.^o"*"*:':)  ^  ^^  points  that  are  paired.  The  line  joining  them 
has  slope  (b^- aj)/[(b2- 02)12  +  (*o~*o)l-  To  find  the  median  slope,  s,  we  compute  the 
median  of  their  inverses,  using  computation  1.  So,  overall,  step  3  takes  O(nlogn)  time. 

Step  4.  Let  s  =  a/[bt2+c]  be  the  slope  found  in  step  3.  A  red  or  black  point,  under  projec- 
tion, parallel  to  a  line  of  slope  s,  onto  the  line  y  =  0,  is  mapped  to  x  =  a-^^t^.  We  need  to 
find  the  (n3+l)-th  (resp  (n4+l)-th)  largest  of  the  projected  red  (resp.  black)  points;  we  do 
this  in  O(nlogn)  time  (computation  1).  Finally,  we  need  to  order  the  two  points  we  found: 
this  requires  determining  the  sign  of  an  expression  ct2+d,  which  takes  a  further  0{n)  time 
(computation  2). 

Step  5.  The  crossing  points  have  the  form  at^  +  b.  We  need  to  find  the  median  of  these 
points;  we  use  computation  1,  taking  time  O(nlogR). 

Step  6.  We  need  to  compute  the  (/i4+l)-th  (resp.  (n3+l)-th  )  largest  of  the  slopes  of  the  set 
of  lines  joining  the  black  (resp.  red)  points  to  p .  Each  of  these  lines  has  a  slope  of  the  form 
a/[btl-^c].  So  we  apply  computation  1  to  the  inverses  of  the  slopes,  as  in  step  3.  This  deter- 
mines the  line  K.  To  finish  the  step  we  proceed  as  follows.  We  project  the  red  (resp.  black) 
points  parallel  to  K,  onto  the  line  ^  =  0;  then  we  sort  the  projected  points  and  the  point  at 
which  K  crosses  >  =  0,  using  computation  1.  We  can  then  determine  those  red  (resp.  black) 
points  to  the  left  and  right  of  K.  So  step  6  takes  a  total  of  O(nlogn)  time. 

Step  7  is  die  same  as  in  the  2-dimensional  algorithm. 

Thus,  Ae  3-dimensional  algorithm  runs  in  O(nlogn)  time. 

In  d  dimensions,  instead  of  applying  the  algorithm  twice  (as  in  the  3  dimensions),  we 
apply  the  algorithm  2**"^  times:  once  for  each  pair  of  regions  with  opposite  signatures  (each 
such  pair  of  regions  is  bisected  by  a  unique  level  d  plane).  Also,  instead  of  having  to  sort 
terms  of  the  form  arJ-t-6,  we  will  have  to  sort  terms  of  the  form  a2tl+  ■  ■  •  +  a^-i'd-i "•"'»<<• 
We  sort  these  terms  using  the  algorithm  for  5(2;),  given  above,  where  we  make  the  following 
change.  For  d>3,  we  note  that  the  problem  Sil;d)  is  effectively  a  problem  of  lower  dimen- 
sion; so  we  solve  this  problem  recursively  (using  the  (<f-l)-dimensional  version  of  the  algo- 
rithm described  here).    It  is  easy  to  deduce  that  the  running  time  for  the  partitioning 
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algorithm,  in  d  dimensions,  d>2,  will  be  0(nlog^('')n),  where /(<f)  =  2"*"^- J. 

Remark  1:  We  can  obtain  similar  partitioning  algorithms  for  a  variant  of  the  partition 
described  in  [YY].  In  [YY]  the  partition  in  J-dimensions  consists  of  a  bisecting  plane  P, 
orthogonal  to  the  axis  x^,  together  with  recursive  partitions  of  the  projection  (in  the  subspace 
P)  of  the  point  sets  in  the  two  halfspaces  defined  by  P.  The  recursive  partitions  are  con- 
strained by  requiring  the  'center'  to  be  equal  in  the  two  partitions  (see  [YY]  for  a  definition 
of  the  center).  However,  it  is  not  required  that  the  same  axes  be  used  in  the  two  recursive 
partitions;  we  impose  this  additional  requirement.  This  reduces  the  number  of  planes  form- 
ing the  partition  to  2**"^. 

To  obtain  algorithms,  similar  to  those  given  above,  we  need  to  reprove  the  existence  of 
the  [YY]  partition,  in  a  style  akin  to  that  of  section  3.  Essentially,  the  [YY]  construction 
requires  the  equating  of  more  and  more  centers,  rather  that  of  angles.  To  show  that  this  is 
possible  we  need  a  lemma  similar  to  lemma  5.  We  state  the  result  below.  We  do  not  state 
or  prove  the  analogue  of  lemma  6;  this  is  relatively  straightforward  and  is  left  to  the 
interested  reader. 

Lemma  7:  Lttfi(x,„...,x{),  f2{xf^,...,x^  be  continuous,  strict,  unbounded  functions,  with/^ 
and  /2  increasing  functions  of  jc^,  ...,  xi,  fi  increasing  in  x^,  and  /j  decreasing  in  x^.  Let 
g(aci,...^2)  ^  ^®  function  such  that  [fi~f2){x^,--.,X2yg)  =  0.  Then/|(x4,...,X2.f)  »s  a  con- 
tinuous, strict,  unbounded  function  of  each  of  Xi^,  ...,  ^2>  ^'^^  >^  '^  ^^  increasing  function  of 
each  of  these  variables. 

We  use  algorithms  very  similar  to  those  given  above;  the  main  difference  is  that  we  pair 
regions  differently  (this  corresponds  to  the  different  pairing  in  the  [YY]  partition).  Also, 
when  seeking  an  angle  rj,  rather  than  seeking  to  make  tj.-^  and  tj^^  equal,  we  will  be  trying 
to  make  two  centers  equal.  But  these  are  differences  of  detail;  the  structure  of  the  algorithm 
is  identical.  The  algorithms  we  obtain  have  the  same  running  time  in  dimension  d  as  those 
given  above,  for  all  d.  The  specific  details  are  left  to  the  interested  reader. 

Remark  2:  An  interesting  question  is  whether  there  are  linear  time  algorithms  for  computing 
the  parallel  planes  partition  in  d  dimensions,  analogous  to  the  linear  time  algorithm  in  2 
dimensions. 

Remark  3:  An  important  feature  of  our  proof  of  separability  is  that  it  has  shown  the  parti- 
tion is  separating  for  the  point  set  at  hand,  rather  than  for  a  continuous  analogue  of  the  point 
set,  as  was  done  in  [C1,YY,DE1],  for  example.  This  leads  to  an  algorithm  naturally, 
whereas  in  the  other  papers  just  mentioned  no  efficient  algorithms  are  given  (essentially, 
those  papers  suggest  using  exhaustive  search).    Our  advantage  arises  because  we  do  not  face 
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the  difficulty  of  having  to  translate  a  proof  about  a  continuous  structure  into  an  algorithm  for 
a  discrete  structure. 
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